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Abstract 

In this report, we present the observation of the exclusive decays B° — ► r] c K° and B + — * r] c K + , 
and the measurement of the related branching fractions. Using a sample of 22.7xl0 6 T(4S) — > BB 
decays collected with the BABAR detector at the SLAC PEP-II asymmetric B Factory during 1999- 
2000, we have observed statistically significant signals in the r\ c — * K^K^ir^ and K + K~ir° channels 
and set upper limits in the ij c — » K + K~K + K~ channels. All the results presented are preliminary. 

We have measured 

B(B+ -> r] c K+) = (1.50 ± 0.19 ± 0.15 ± 0.46)xl0" 3 
B{B° r) c K°) = (1.06 ± 0.28 ± 0.11 ± 0.33)xl0- 3 

where the first error is statistical, the second systematic and the last due to the uncertainty on the 
world average r\ c — ► KKtt branching fraction. 
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1 Introduction 

We present the measurement of the branching fractions of the exclusive decays^ B° — > rj c K° and 
B+ rj c K+, with 7] c decaying into K%K ± ir f: , K+K'tt , and K + K~ K + K~ {K° 3 vr+vr- and 
7T° — > 77). The 7] c is a cc meson with I G {J PC ) = 0+((T + ). The decay B° -> t? c K° proceeds 
through the same b — ► cc s color-suppressed quark diagram as the "golden" mode, B° — ► J/tpK°, 
used to measure the CP-violating parameter sin2/3 with negligible theoretical uncertainty jl]]. Up 
to now, experimental information on B decays into ry c has been sparse Q. 
The ratio of the decay rates for the exclusive charmonium decays 

R K = T(B^r, c K)/T{B^J/^K) (1) 

has been calculated with different dynamical assumptions 0]— [§j including factorization^]. The 
ratio is used since one expects that the corrections to the heavy quark limit, due to the relatively 
light s-quark, are likely to cancel. This leads to the following predictions for Rk'- 1-6 ± 0.2 Q, 
1.64 ± 0.55 §, 1.8 ~ 2.3 §, 0.94 ± 0.25 0, 1.0 ~ 1.3 §. 



2 The BABAR detector and dataset 

The data used in this analysis are obtained with the BABAR detector at the PEP-II asymmetric 
e + e~ storage ring. The BABAR detector is described elsewhere [ 10 ] . The 1.5 T superconducting 



solenoidal magnet, whose cylindrical volume is ~ 1.4m in radius and ~ 3m long, contains a charged- 
particle tracking system, a Cherenkov detector (DIRC) dedicated to charged particle identification 
and an electromagnetic calorimeter. The segmented iron flux return, including endcaps, provides 



identification of muons and K®. In addition, the end of the cylindrical volume in the e~ direction 
is instrumented with an electromagnetic calorimeter. The tracking system consists of a 5-layer 
double-sided silicon vertex tracker and a 40-layer drift chamber filled with a gas mixture of helium 
and isobutane. The calorimeter consists of 6580 CsI(Tl) crystals. The flux return is instrumented 
with resistive plate chambers. 

We have used data corresponding to 20.7 fb _1 of integrated luminosity collected at the T(4S) 
resonance ( "on-resonance" ) , and 2.1 fb _1 recorded ("off-resonance") about 40MeV lower in energy 
in the T(45) rest frame ("CM"), between October 1999 and October 2000. The asymmetric 
collisions produce a boost in the e~ direction, with ^7 = 0.55 in on-resonance running. 

1 Throughout this paper, whenever a mode is given, the charge conjugate (cc.) is also implied. 

2 We note that we have measured a departure from the factorization hypothesis Q in another b — > cc s color- 
suppressed mode, B — » J/tp K* , wherein we have made a polarization measurement, more sensitive than a measure- 
ment of Rk to the existence of a factorization-violating term. 
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3 Analysis method 



A blind analysis is performed in which all selections are chosen to maximize N$ / \J N$ + Nb using 
simulated or off-resonance data, or sidebands in on-resonance data. Ns(Nb) is the number of 
expected signal (background) events after all selection criteria have been applied. 

Event selection designed to enhance the number of B decays requires four or more charged 
tracks, the sum of all charged and neutral energies to be above 2 GeV, the sum of all the charged 
momenta to be above 1 GeV/c, and the normalized second Fox- Wolfram moment |jll[ to be less than 
0.6. In addition, at least one neutral or charged kaon candidate is required to have a momentum 
in the T(4S) rest frame consistent with the two-body decay B — > r\ c K. 

B° or B^ candidates are formed from an r\ c candidate and a "fast" kaon, either a charged kaon 
or a Kg —* 7r + 7r~. The r\ c candidates correspond to three different topologies: two charged tracks 
with either Kg — > tt + tt~ or ir° — > 77, or four charged tracks. The B decay vertex is calculated 
using the charged r\ c daughters, and the fast kaon if charged. 

We require that all charged tracks be within 0.35 < 9 < 2.54 to obtain well-reconstructed 
tracks, where 6 is the polar angle with respect to the e~ direction. An important requirement of 
our analysis is that charged kaon candidates from the B — > rj c K decay are identified by the DIRC 
and/or by measurements of ionization energy loss dE/dx in the drift chamber and silicon tracker. 
The momentum of each kaon from r\ c decay is required to be greater than 250MeV/c. 

The Kg particles can arise from r\ c or B decays. In the following, the number in parentheses 
corresponds to the latter Kg. The Kg — > it + -k~ candidates are required to have a reconstructed 
invariant mass within 12.5 (10) MeV/c 2 of the nominal, i.e. world average, value [p^| . Furthermore, 
the cosine of the opening angle between the flight direction and the momentum vector of the Kg 
candidate is required to be greater than 0.990 (0.9995), and the flight distance from the B vertex 
greater than 2 (3) times its error. 

The 7T° — > 77 candidates are formed from pairs of photons detected in the calorimeter with a 
reconstructed invariant mass within 15 MeV/c 2 of the nominal value. We require that the cosine of 
the decay angle in the 7r° rest frame be less than 0.82 to avoid accidental combinations involving 
very soft photons. In addition, the electromagnetic showers are required to have moments of the 
lateral energy deposition |ll|] between 0.01 and 0.55. The lower energy photon has a minimum 
energy of 130 MeV while the minimum value for the higher energy photon is 270 MeV. 

A Fisher discriminant is used to suppress continuum backgrounds. The Fisher variable is 
defined as a linear combination of eighteen variables, including the energies between each of nine 
cones relative to the r\ c direction in the CM |jl4j| . The most important variables are the normalized 
second Fox- Wolfram moment and the event thrust, constructed with all charged tracks and neutral 
clusters in the event. The Fisher discriminant is trained on signal, uu, dd, ss, and cc simulated 
events, and tested on off-resonance data. The requirements on the Fisher discriminant depends on 
the decay mode. 

The charmonium mass region is defined by 2.74 < rax < 3.22 GeV/c 2 . After all selection 
criteria, the weighted double- Gaussian mass resolutions are 10, 12, and 26 MeV/c 2 for the KgK^ir^, 
K + K~ K + K~ and K + K~ir° channels, respectively, as obtained from a simulation. The rj c signal 
region varies between ±55 MeV/ c 2 and ±70 MeV/ c 2 relative to the nominal r\ c mass (2979.8 MeV/ c 2 ), 
depending on the r\ c decay mode. 

The total energy of the e + e~ system in the T(45*) CM and laboratory frames are denoted 
by and E Q , respectively. In the e + e~ laboratory frame, the candidate energy is defined as 
Eb = (s/2 + p G • pb)/E q , where p D and ps are the momentum vectors of the e + e~ system and 
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the B candidate, respectively pi . The analysis region is defined by a rectangular area in the AE- 
tties plane where AE is the difference between the energy of the B candidate in the CM frame 



and \/s/2, and tties is the beam-energy substituted mass, d E 2 B —p 2 B . For events with multiple 
candidates, the one with the smallest \AE\ is retained; this choice affects only a small fraction of 
events, from 3.4% to 12.4% in the analysis region. 

The limits of the analysis region are defined by 5.1<rriES <5.29 GeV/c 2 and \ AE\ <0.25 GeV. 
According to the full detector simulation based on GEANT3 [13[, depending on the r\ c decay mode, 
the signal is Gaussian-distributed in AE with a mean near zero and a resolution between 15 and 
30MeV, while it is Gaussian-distributed in tties with a mean near the B mass and a resolution 
around 2.5 MeV/c 2 . The AE resolution depends on the rj c decay mode, best for K^K^tt^ and worst 
for K + K~7T°. Note that the AE distribution in data is not centered at zero but rather at about 
— lOMeV; the window is shifted accordingly, leading to a contribution to the overall systematic 
error. The shifted- 1 AE\,\m E s - m B ominal \ si gnal region is <30 MeV, <7 MeV/c 2 for the tightest 
(i^°K ± vr =F ) case and <70 MeV, <9 MeV/c 2 for the loosest (K+R-tt ) case. 



4 Observation of exclusive r\ c signals 

Figure p] displays the mass distribution of the charmonium system in the (AI^mEs) signal region 
for the KgK^n^ channel, using B + candidates after subtraction of the combinatorial background. 
We see clear r\ c and J/ip peaks where we have indicated the rj c mass selection excluding the J/ip 
region. The representative curves are fits of three contributions: flat background, J/ip peak, and rj c 
peak with two different widths. The J/ip peak is represented by a Gaussian with mean constrained 
at the nominal J/ip mass and a 12 MeV/c 2 resolution. The rj c mass peak is represented by a 
Breit-Wigner distribution convoluted with the same Gaussian. The mean of the Breit-Wigner 
distribution is fixed at the nominal rj c mass and the width^, either to the world average or to the 
CLEO measurement [16]. Since we cannot yet distinguish among the various measurements, we 
have used for the r\ c width the average value 16.7 ± 6.0 MeV/ c 2 ; the efficiency depends on the width 
central value and the systematic error on its error. 

In Figure ^ we display the analysis region for the B + and B° (rj c — > K^K^tt^) channels as 
examples. Clear accumulations in the (AE,m,Es) signal region are apparent. Figures || and |4| display 
projections of the analysis region for the different r\ c channels. The combinatorial background shape 
is parametrized by a linear function in AE and a threshold function |l^] in m^s with a fixed endpoint 
given by the average beam energy. 

In addition to the combinatorial background, a background that peaks in the (AE,m,Es) signal 
region can arise from cross- feed from other r] c decay modes, from partial reconstruction and/or 
incorrect particle identification, or from B decays into the same detected particles without an 
intermediate r] c decay (exact matches). After study, the first two sources are found to be negligible. 
A quantitative evaluation of the exact matches for each mode is made using data by studying the 
r) c mass sidebands for events in the (A£,?tjes) signal region, after subtracting the combinatorial 
background as a function of mass. The peaking background is consistent with zero for all modes 
except possibly for the K^K^ir^ mode, see the flat background in Figure |l[ 

The raw yield and expected backgrounds in the (A.E,mEs) signal region, and the probability 
that the background fluctuates to the observed yield are given in Table [l]. In order to ensure the 

3 The world average width is 13. 2±|1 MeV/c 2 while more recent results give 27.0 ± 5.8 ± 1.4MeV/c 2 (lfi), 11.0 ± 
8.1 ±4.1 MeV/ c 2 and 21.1^MeV/c 2 ft| 
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Figure 1: K°K mass for B + candidates in the (AE, mEs) signal region after subtraction of the combinatorial 
background. The fits are described in the text. The remaining flat background is that due to the peaking background; 
see text. The "thick" curve corresponds to the CLEO 7) c width while the "thin" curve corresponds to the world 
average. 
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Figure 2: AE vs. m-ES for candidate B + — > r) c K + events (on left) and B° — > rj c Ks events (on right), with 
r/ c — ► K° s K ± ir 1L . The (AE,m-Es) signal region is indicated. All selection criteria have been applied except for the 
signal region requirements. 
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Figure 3: The AE distribution relative to its mean in the rriES signal band for combined B + — > r) c K + and 
B° — > T) c Ks candidates with r\ c — > K%K n^ 1 , fitted to a double Gaussian with common mean on top of a linear 
background. The weighted average resolution is 16.3 MeV/c 2 . The narrower Gaussian represents 71% of the area of 
the double Gaussian; its resolution is 5.9 MeV/c 2 . All selection criteria have been applied except that for AE. 



statistical independence of the signal and background measurements, the combinatorial background 
is estimated here by the extrapolation into the toes signal band of the threshold function fitted in 
the AE signal band below the toes signal band (toes < 5.27GeV/c 2 ). Because of the low statistics 
in the B° channels, the shape parameter of the background function is fixed to that fitted in the 
corresponding B + channel. 

Table 1: Raw yield, extrapolated combinatorial (see text) and peaking backgrounds in the (AE,m-Es) signal region, 
and Poisson probability that the combined background fluctuates to the number of events found in the signal region 
(called "Prob/!„ ct "). Due to the limited data sample, the fitted combinatorial background estimate for B° (;q c — * 
K + K~ K + K~) comes from the AE sidebands. 



mode 


Yield 


Fitted 
combinatorial 
background 


Peaking 
background 


Prob/ itlct 


B + ( Vc K2K=ir + ) 


72 


6.08 ± 1.39 


6.12 ± 2.61 


2 x 10~ lb 


B + (j) c -> K+K-tt u ) 


25 


2.92 ± 0.92 


0.58 ± 0.58 


3 x 10" i& 


B+ (r) c -»■ K+K'K+K-) 


17 


7.41 ± 1.78 


1.72 ± 2.75 


2 x 10~ a 


B" ( V c K v s K=n + ) 


19 


1.18 ± 0.38 


1.48 ± 1.08 


3 x 10 _1J 


B {> { Vc -> K+K-ir u ) 


8 


1.73 ± 0.38 





4 x 10~ 4 


B u (r? c -f K+K-R+K-) 


1 


1.01 ± 0.25 
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Figure 4: mes distributions in the AE signal band for candidate r\ c decays: K%K ±r K^ (top left), K + K~ty° (top 
right), and K + K~ K + K~ (bottom). All plots use combined B + and B° data and are fitted, using a binned likelihood 
method, to the sum of a Gaussian signal, with an average resolution of 2.9 MeV/c 2 , and a threshold background 
function Jl9|. All selection criteria have been applied except that rnES- 



5 Branching fraction determination 

The measured B + or B° branching fraction (J3) is given by 



B ~N BS Xe> (2) 
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where N y i e id is the net yield in the (AE ,m^s) signal region, extracted from fits to the tbes distribu- 
tions in the AE signal region (Figure ||), and corrected for the peaking background contributions 
listed in Table [lj; e is the signal efficiency determined by applying the same analysis chain to signal 
Monte Carlo (MC) samples and correcting for data-MC differences; and N BB is the number of 
produced BB pairs, (22.73± 0.36) xlO 6 , determined by a comparison of the rate of multihadron 
events taken on-resonance to that off-resonance. 

5.1 Determination of signal efficiency 

The efficiency for reconstructing B — > r/ c K candidates for each r\ c decay mode is given by the fraction 
of generated signal events that are reconstructed in the appropriate mode. We have compared 
simulations with real data, using for example t + t~ and D*^ control samples. There are small 
differences in reconstruction efficiency for charged particles, K® and 7r° mesons, vertexing efficiency, 
resolution and absolute scale of charged particle momentum and photon energies, and charged 
kaon identification and pion misidentification probabilities. These effects have been measured and 
corrected. The resulting efficiencies are given in the first line of Table |2[ 

Table 2: Relative systematic errors on efficiency. All values are expressed in percentage relative to the efficiency, 
which is given in the first line as a fraction. The last line gives the total relative systematic error obtained as a 
sum in quadrature of the individual contributions. The 1.6% error from the determination of the number of BB 
events, common to all modes, is not listed but is included in the total as is the statistical error on the efficiency 
determination. 



t) c decay 


K+K-K+R- 


B u - IMS! 


K+K-n° 


B+ Vc K+ 
K+K-K+K- 


B+ -» r lc K+ 


B+ -> Vc K+ 


Efficiency 


0.111 


0.148 


0.0733 


0.117 


0.145 


0.0635 


Rel. stat. err. 


0.004 


0.003 


0.0027 


0.003 


0.003 


0.0017 


Tracking eff. 


9.5 


7.7 


5.8 


7.8 


6.2 


4.4 


K" eff. and cuts 


5.9 


12.2 


5.3 




6.8 




7 eff. and ty" cuts 






3.5 






3.5 


Vertexing eff. 


1.3 


1.0 


1.0 


1.3 


1.0 


1.0 


Kaon ident. eff. 


10.5 


2.9 


5.6 


14.1 


6.5 


9.2 


Fisher cut eff. 


2.3 


2.2 


4.0 


1.1 


2.9 


4.0 


77 c width uncert. 


2.9 


3.2 


3.1 


2.9 


3.2 


3.1 


AE centroid shift 


0.47 


3.3 


0.86 


0.51 


2.7 


0.24 


AE resolution 


3.6 


3.8 


4.4 


3.4 


5.2 


3.4 


E 


16.2 


16.1 


12.3 


16.8 


13.5 


12.4 



5.2 Determination of systematic errors 

We have evaluated the systematic errors on the yield, B counting and efficiency determination. The 
systematic error on the yield comes from a comparison of the combinatorial background estimations 
from the AE side and signal bands while that on B counting comes principally from the uncertainty 
on the efficiency due to small differences between data and simulation. 

Each of the efficiency corrections, as well as our knowledge of the r/ c width, has a corresponding 
systematic uncertainty. In addition, each requirement in the analysis method has been studied to 
evaluate any systematic differences between simulation and data. The dominant systematic errors 



13 



on the signal efficiency are due to kaon identification, tracking efficiency, and K® reconstruction as 
can be seen in Table |. 

5.3 Results 

Our results for the product of the branching fractions for each mode are listed below. We have 
used the nominal values for the K® —* -k + tt~ and ir° — ► 77 branching fractions. The B — > r\ c K 
branching-fraction determinations assume that the branching fraction of the T(4S) into BB is 
100%, with an equal admixture of charged and neutral B final states, and similarly for K° relative 



to K° and K°. 










B(B+ - 


-+ r /c K+)fi(r ?c - 


^K°K 


7T + + C.C.) 


= (52.8± 7.9± 7.3) xl0~ 6 


B{B+ - 


-» r ?c K+)^(r ?c - 




-vr ) 


= (15.5± 3.6± 2.5) xl0~ 6 


B{B+ - 


-» r ?c K+)B(r ?c - 


-> if+if 


-K+K-) 


< 5.6xl0- 6 (90% CL) 


B{B°- 


- r) c K°)B( Vc - 




7T + + C.C.) 


= (36.8± 11.6± 6.0) xl0~ 6 


B(B°- 


-» VcK°)B( Vc - 


> K + K~ 




= (11. 3± 5.1± 2.4)xl0~ 6 


B{B°- 


-» VcK°)B(rj c - 




K+K-) 


< 2.3xl0~ 6 (90% CL) 



The first error is statistical and the second systematic. The central value for B(B + —* r] c K + )B(r/ c — > 
K+K-K+K-) is 3.2 xl0~ 6 , while the two-sided 68% CL varies from 2.6xl0~ 6 to 4.1xl0" 6 . No 
correction is made for any potential contribution to the K + K~ K + K~ channels. 

The channels rj c — ► Jf^lf^Tr^ and r/ c — ► i^ + ir _ 7r° are manifestations of the general decay 
77c — * KKir. From isospin symmetry, the corresponding rates are related by simple Clebsch-Gordon 
coefficients: B{r] c -» K K~Tr + + c.c.) = 2/3 B(rj c -> A^tt) and /3(r? c -» K+K'tt ) = 1/6 Sfe -> 
KKtt). Therefore the ratio of branching fractions, B(rj c — > K + K~tt°) /B(i] c — > K°/f _ 7r + + c.c), 
should be 0.25. Our measurements are consistent with this value for i? + (0.29 db 0.08 ± 0.04) and 
5° (0.31 ± 0.17 ± 0.05). 

We therefore combine our two results, taking into account common systematic errors, to obtain 
the values for the general decay: 

B(B+ -» rj c K + )B{ri c -> AiTvr) = (82.5 ± 10.4 ± 8.3)xl0~ 6 
-> r] c K°)B{r} c -> AAV) = (58.1 ± 15.2 ± 6.3) xl0~ 6 . 

The first error is statistical and the second systematic. We deduce the branching fraction ratio from 
our measurements of the KKtt channel: B{B° r) c K°)/B{B + r? c K + ) = 0.71 ± 0.20 ± 0.08. 
We have not used the 7/ c — ► K + K~ K + K~ results since their statistical weight would be marginal. 

Using the world average for the r\ c — > KKtt branching fraction, 0.055 ± 0.017 our results 
become 

B{B + -» r] c K + ) = (1.50 ± 0.19 ± 0.15 ± 0.46)xl0~ 3 
B(B° -» r) c K°) = (1.06 ± 0.28 ± 0.11 ± 0.33)xl0~ 3 , 

where the last error is due to the r\ c — ► KKtt branching fraction. We have not used the r; c — ► 
If \K'~.K' + .K'~ results since the r\ c branching fraction is not very well known. We compare these 
results to the exclusive branching fractions measur ed by CLEO g: B(B + -► r Jc K+) = (0.69±g:|? ± 
0.08 ± 0.20) x 10~ 3 and B(B° -> i] c K°) = (1.09^;|| ± 0.12 ± 0.31) x 10" 3 The third error is that 
due to the nominal J/0 — > 777c branching fraction. Assuming that the errors due to the nominal 
branching fractions cancel, our results differ by a factor 2.2 ± 0.9 for the B + channel, combining 
statistical and systematic errors in quadrature. The B° channel results are consistent. 



14 



To determine Rk, we have used our measurements |2(J of the branching fractions, B{B + — > 
J/^K + ) = (10.1 ± 0.3 ± 0.5) x 10" 4 and B(B° -> J/^K°) = (8.5 ± 0.5 ± 0.6) x 10" 4 , taking into 
account common systematic errors, to obtain 

R+ = T(B + -> r] c K + )/T(B + -» J/ipK + ) = 1.48 ± 0.19 ± 0.17 ± 0.46 
R0 R = r(B° -> r] c K°)/T(B -» J/V'K ) = 1.24 ± 0.33 ± 0.16 ± 0.38 

where the first error is statistical, the second systematic and the third due to the rj c — ► KKir 
branching fraction. Our results agree with the theoretical predictions listed at the end of Section 
1. 
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